Pairing fluctuations in trapped Fermi gases 
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Fluctuations of the amplitude of the order parameter govern the properties of superconducting 
systems close to the critical transition temperature. In the BCS regime we examine the contribution 
of these pairing fluctuations to the superfluid order parameter for harmonically confined atomic 
Fermi gases. In the limit of small systems we obtain an expression for the space dependence of 
the fluctuations, in good agreement with the results of numerical calculations. In this limit we also 
predict a parity effect, i.e. that pairing fluctuations should show a maximum or a minimum at 
the centre of the trap, depending on the value of the last occupied shell being even or odd. We 
finally propose to detect pairing fluctuations by measuring the density-density correlation function 
as evaluated after a ballistic expansion of the gas. 
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Progress in trapping and cooling of gases of fermionic 
atoms have paved the way to explore a rich variety 
of phenomena. A major experimental effort nowadays is 
directed towards the observation of the superfluid tran- 
sition. The presence of a Feshbach resonance allows to 
study the behaviour of Fermi gases both in the weakly 
and in the strongly interacting regime by varying an ex- 
ternally applied magnetic field. Away from the resonance 
the effective interatomic interaction is well represented 
by a contact pseudopotential gS(r), with g — AirTi^a/m 
and a being the s-wave scattering length. A number 
of papers have explored the properties of a Fermi sys- 
tem across the resonance with particular attention to the 
BCS-BEC crossover and the strongly-interacting limit 
kpa — > ±oo with kp = (37m) 1 / 3 and n being 
the atomic density. In the regime kpa —> + bosonic 
molecules have been already observed Q and they have 
been reported to undergo condensation 0; current ex- 
periments start to explore the a < region It has 
been suggested that by adiabatically tuning the scatter- 
ing length from positive to negative it could be possible 
to reach a BCS-like regime Here we concentrate on 
such weak-coupling regime kpa — > - Q. 

Mean-field BCS theory predicts the appearence of a 
non- vanishing order parameter Abcs (T) below a critical 
temperature T c which depends strongly both on the in- 
teraction strength and on the geometry of the system || . 
The scale of variation of the order parameter is the coher- 
ence length £ = h 2 kp j mir Abcs (0) 3> n~ 4 / 3 , m being the 
atomic mass, to be compared with the size R of the cloud 
under confinement (which we suppose harmonic). In the 
case £ <C R of very large traps mean-field theory is very 
accurate in describing the physical properties of the con- 
densate. The experiments with atomic gases, however, 
are performed on systems of adjustable size ranging from 
clouds of about 10 4 — 10 7 atoms in harmonic traps to 
droplets of few tens of atoms in each lattice site of an op- 



tical lattice. One expects that by lowering the number of 
atoms in the trap (or equivalently when £ > R) fluctua- 
tions of the order parameter can become important close 
to the critical temperature. This is the regime we dis- 
cuss in this work. Above T c fluctuations signal the onset 
of the Cooper pair instability which will eventually drive 
the system to the superfluid state. In superconductors 
the role of fluctuations both on the thermodynamics and 
on the transport properties has been studied in great de- 
tails both theoretically and experimentally and we refer 
to Ref. Q for a comprehensive review of the field. 

The aim of this Letter is to formulate the theory of 
pairing fluctuations in the case of trapped Fermi gases. 
This problem resembles that studied in zero-dimensional 
superconductors (i.e. of radius smaller that the coherence 
length) P3- 

There are however important differences 
that we are going to highlight. First of all pairing fluctu- 
ations are inhomogeneous due to the presence of the trap; 
we will show below that they are enhanced in the center of 
the trap. Another important issue is how to detect fluctu- 
ations. In conventional superconductors they dominate 
the behaviour of most thermodynamical and transport 
properties (e.g. paraconductivity, density of states, spe- 
cific heat, . . . ) close to T c 9]. In trapped atoms, however, 
those quantities are very hard or impossible to measure. 
As a clear signature for the detection of pairing fluctua- 
tions we propose to analyze density-density correlations 
which were recently suggested by Altman et al. as a 
probe of many-body states in ultra-cold atoms. 

We assume that the system is composed of fermions 
in two internal states a —1,1 with equal numbers and 
attractive intercomponent interactions, which we shall 
describe by the contact pseudopotential gS(r). Corre- 
spondingly, we shall consider only pairing fluctuations in 
s-wave. This is the main collisional channel at ultralow 
temperatures, while s-wave interactions among atoms be- 
longing to the same internal state are forbidden by the 
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Pauli principle. The second-quantized Hamiltonian reads 

<r=U \ m j 

+ g [ d 3 r^(r)^(r)^(r)V T (r), (1) 



where ip a are the fermionic field operators, /x is the chem- 
ical potential, V(r) the confining potential and g < 0. 
Close to T c the most important contribution beyond 
mean field comes from pair fluctuations, we therefore 
treat the density correlations at the Hartree level and 
add the Hartree fields W a {r) = gn_ CT (r) to the single- 
particle part of the Hamiltonian. 

By introducing the auxiliary field A through a 
Hubbard- Stratonovich transformation |l^ the fermion 
fields can be integrated over. Close to the criti- 
cal temperature the partition function becomes Z = 
Z J d[A*]d[A]e- s, [ A *< A ]/ ?i , where Z is the partition 
function of the gas in the Hartree approximation and 
the effective action to quartic order is 

S = -hpj dndra A*(n)A(n,r 2 )A(r 2 ) (2) 
+ T / II ^A*(r 1 )A*(r 3 )i?({r l })A(r 2 )A(r 4 ), 

•* j=l,4 

where [3 = 1/fesT and we have ignored the dependence 
of A on the imaginary time since we are interested in the 
classical fluctuations above T c . The kernel j4(ri,r 2 ) — 
g -1 5(ri — r 2 ) + h~ 1 K(ri, r 2 ) is defined through the two- 
particle propagator 

1 ^ tanh(/3&/2) + tanh(/3^-/2) 



K(r u r 2 ) = 



& + & 



x ^(ri)^-(ri)0?(r 2 )^*(r 2 ) 



(3) 



expressed in the Hartree- Fock basis {<pi(r)} with ener- 
gies {ti} (£i — e,- — fi). The kernel K requires a proper 
regularization as implied by the use of the contact pseu- 
dopotential |l3j . The kernel of the quartic term is given 
by 

B({Ti}) = n f H ^60(1,4)00(1,2)60(3,2)00(3,4). 

J i=l,3 

(4) 

where Go (a, b) is the single-particle Green's function and 
the variable label a — 1, .., 4 stands for (r a , r a ). 
Pairing fluctuations - In order to evaluate the pair fluc- 
tuations (|A(r)| 2 ) = Z' 1 J d[A*]d[A] \ A\ 2 e -^ A '' A ^ h it 
is convenient to decompose A(r) in normal modes 

A(r)=£ X „(r)A„, (5) 

where Xi/( r )'s are the eigenvectors of the kernel A 0] 
Sri A(t, ri)xiz(ri) = a„x^(r). (6) 



For temperatures larger than T c the dominant contribu- 
tion to the fluctuations comes from the quadratic part 
of the effective action and in particular, from the lowest 
eigenvalue ao (and the associated eigenvector Xo( r )) 



(|A(r)| 2 



a 



Xo(r)| ; 



(7) 



For a generic system under confinement both ao and 
Xo( r ) have to be calculated numerically by diagonaliza- 
tion of the kernel A. 

For an isotropic harmonic trap V(r) = muj 2 r 2 /2 we 
proceed now to derive analytical results in the small- 
system limit, i.e. when Rtf "C £, where Rtf = 
(2ep /moj 2 ) 1 / 2 is the Thomas Fermi radius of the gas 
and £ the bulk BCS coherence length estimated at the 
center of the cloud. Such condition is equivalent to 
£ sT buik ^ ^ where T c bulk is the bulk prediction for 
the critical temperature of the BCS transition evaluated 
using the central density of the cloud. In this limit the 
pairing energy is much smaller than the trap level spacing 
and the Cooper pairs are only formed between particles 
with quantum numbers (n,l,m) and (n,l,—m) residing 
in the same harmonic oscillator shell ( "intrashell pairing 
regime" ) [|| . Since for an isotropic trap the lowest eigen- 
mode is spherically symmetric the eigenvalue equation 
© becomes 0| 



[1 - \g\ao(T)]Xo(r) = J dr' r' 2 gK(r, r')xo(r'), (8) 

where the kernel describing spherically symmetric solu- 
tions is 



K(r,r>) = ~ ]T(2; + 1 



tanh(/?£„i)+t,anh(^„<//2) 



n,n' ,1 

x R nl {r)R n ,i{r)R nl {r')R nn {r') 



(9) 



In the small-system limit the effect of the Hartree field is 
negligible and the R n i(r)'s are approximately the radial 
eigenf unctions of the isotropic harmonic oscillator. Cor- 
respondingly £ n i ~ (n + 3/2)hto — (i, where /1 ~ cf = 
(rip + 3/2)Klu is the chemical potential. We further as- 
sume that pairing strictly only occurs within the same os- 
cillator shell, and that it mostly occurs about the Fermi 
surface. In particular, we simplify Eq. JHJl by setting 
/ dr' r' 2 R n i(r')R n >i(r')xo(r') ^ K,n'Xn 0, by assum- 
ing a weak dependence on n for the relevant shells around 
the Fermi level, and by taking angular averages. The 
above assumptions lead to the final expression for the 
normalized eigenfunction Xa(r) and the corresponding 
eigenvalue ao which determine the spatial and temper- 
ature dependence of the pairing fluctuations in Eq. J7J , 



Xo(r) = 



15 {ira scY /2 y^ (2/ + l) 
2 (2n F + 3) 5 / 2 2f 4tt 



R nF i(r)\ 2 (10) 
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a (T) = 



4 (3 c Rtf 
157T 3 a* 
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(11) 



where T c has been computed in Ref. 8], a osc — 
(Ti/moj) 1 / 2 is the harmonic oscillator length and we have 
kept only the leading temperature dependence in ao- 
An explicit expression for Xo(r) can be obtained by us- 
ing that J2i( 21 + l )\ R nA r )\ 2 l^ = dp{r)/dn F , with 
p(r) being the single-component density profile. In the 
Thomas- Fermi approximation for p(r) this yields Xo( r ) = 
(15/8ttR^ f ) 1 / 2 (l - r 2 /i?T F ) 1/2 > showing that the gas is 
mostly susceptibile to pairing fluctuations in the central 
part of the trap. 




FIG. 1: Lowest eigenvector Xo(f) (in units of aosJ 2 ) as a func- 
tion of the radial coordinate r (in units of a oac ), as obtained 
from numerical diagonalization (solid line), compared with 
the analytical solution HUH (dot-dashed line) and with the 
Thomas- Fermi approximation (dashed line). The inset shows 
(|A(r = 0)| 2 ) 1//2 (in units of fcsT c ) as a function of tempera- 
ture (in units of T c ) as obtained from the quartic (solid line) 
and the quadratic approximation (dashed line) of the effective 
action, and the mean-field solution (dotted line). 



The spatial behaviour Xo( r ) of the pairing fluctua- 
tions is illustrated in Fig. ^ f° r a system with n F — 
10 filled shells (~ 500 particles) and coupling constant 
\g\ = 0.7 huja^ sc , as obtained from the analytical expres- 
sion (|10|l as well as from the numerical results obtained 
by diagonalization of the regularized kernel A with a pro- 
cedure analogous to that described in Ref. 01 • F° r our 
choice of parameters ksTc/huj = 0.038 <C 1, and the 
system is well within the intrashell regime (T c ~ InK). 
We find that the analytical expression well approximates 
the numerical solution, while the characteristic oscillat- 
ing behavior due to the discrete level structure is lost in 
the Thomas-Fermi approximation. From the full quan- 
tum shell-structure calculations we also predict a parity 
effect, i.e. that pairing fluctuations should show a maxi- 
mum or a minimum at r = 0, depending if n F is even or 
odd. 

At T = T c the lowest eigenvalue ao vanishes, indicat- 
ing the onset of the superfluid phase. The temperature 




FIG. 2: Lowest eigenvalue ao(T) in units of (fesT c )~ 1 a~c 
as a function of temperature (in units of T c ) as calculated 
numerically (x's) and analytically via Eq. ill 1 1 (solid line). 
The inset shows the lowest eigenvalues a v (in the same units 
as ao) as functions of the index u, evaluated at T ~ T c . 



behaviour of ao is illustrated in Fig. El where we compare 
the approximate expression Ijlll) with the exact numerical 
solutions for the same choice of parameters as in Fig. ^ 
The figure shows an excellent agreement between the two 
calculations. 

Since ao(T c ) = 0, according to Eq. Q the pairing fluc- 
tuations diverge at T c (dashed line in the inset of Fig. ^| . 
This is an artifact due to the quadratic approximation 
used for the effective action, and can be cured by taking 
into account the quartic term [l(ij |. which we estimate 
from Eq. (J2J within a single-mode approximation. This 
is well justified by noticing that the lowest eigenvalue is 
significantly smaller than the other ones (inset in Fig.|2J, 
and thus the v = mode strongly dominates the fluctu- 
ations in the intrashell regime. The effective action then 
reduces to the simple form 



S e s — fofic 



a |A,=o| 2 + -B (T C )|A„ =0 | 4 



(12) 



where B (T C ) ~ 64/3j?/6757r 6 n_Fao SC and we have only 
retained the contribution to the Green's functions in 
Eq. J2J coming from the levels at the Fermi sur- 
face, consistently with the intrashell pairing approxi- 
mations. By introducing the dimensionless quantity 
A„ = o/(«osc ksTc), Seff only depends on np and T/T c , 
since RTF/a osc = (2n F + 3)^ 2 . The value of riF controls 
the magnitude of the critical region around T c , which we 
estimate to be 5T/T C - 2/(n F + l)(n F + 2) from the 
Ginzburg criterion. 

Using the effective action Eq. I|12fl the temperature be- 
haviour of the pairing fluctuations can be estimated an- 
alytically, and is illustrated in the inset of Fig. Q] (solid 
line). The figure shows that the divergence predicted by 
the quadratic term alone is indeed eliminated and the 
behavior is now smooth across T c , tending to the mean- 
field result (dotted line) below T c . We should remark 
that Eq. I|12|) is only accurate close enough to T c , while 
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well below T c the full action should be withheld in the 
averaging. 

We also mention that in the large-system limit £ <C 
Rt f pairing fluctuations can be described along the lines 
of Ref. 01- We do not pursue the issue further in this 
work since in this limit they are expected to be small. 
Density correlations - The presence of Cooper pairs im- 
plies non-zero density-density correlations among parti- 
cles with opposite spins located at q, q. Altman et 
al. 0| suggested that such "anomalous" correlations 
could be measured probing the density-density correla- 
tion function after the trap has been released and the gas 
has undergone a ballistic expansion: (n a (r)n_ a (r')) ~ 
(n cr (q r )n_ (T (q r /)), where q r = (m/M)r, m is the mass of 
the fermions and t is the time lapse after the trap has 
been opened. At mean-field level it can be shown that 
(5n a (r)5n- a (-r)) oc \A B cs\ 2 , with Sn a (r) = n a (r) - 
(n CT (r)), and hence, one could directly access the gap 
field amplitude. Above T c the same density-density 
correlation function is also sensitive to pairing fluctua- 
tions. Proceeding as for obtaining Eq. J2J) for a uni- 
form system we get (<5n (T (r)(5n_ (T (— r)} = (((A/2_E<j r ) 2 [l — 
2f(E qr )] 2 ))s afs , where ((-))s e// means the average over 
the distribution of A weighted by the effective action 
S e ff, f(x) — {exp[/3(x — u,)} + is the Fermi distri- 
bution and E q = {fi A q A /Am 2 + A 2 ) 1 / 2 is the usual BCS 
quasi-particle dispersion relation. 

In small harmonically trapped systems anomalous cor- 
relations occur between particles with opposite spins and 
quantum numbers (n,l,m) and (n,l,—m). In the in- 
trashell limit analogous approximations to those used in 
order to obtain Eqs. (fTUfl and ijTT|) yield 

(5n a (v)Sn^(-v)) = I {^^j «[F(A)] 2 )) Sc// , 

(13) 

where F(A) = 1 - 2/(|A|) + 2(\A\/Tilu) log( 7 n F ), with 
7 = e c ~ 1.78, and C Euler's constant. In this limit 
the spatial dependence of the anomalous correlations is 
proportional to |xo( r )| 2 j which is given by Eq. i|10|) and 
shown in Fig. ^ Also for the small system the mean- 
field result is non-zero only below T c , where we find 
(5n ff (r)<5rc_ CT (-r)) oc |A BCS (r)| 2 . 

The effect of pairing fluctuations is not confined to the 
density-density correlations described in this work. As 
another example, we expect a soft gap in the laser in- 
duced tunneling (discussed in Ref. [l||) above the critical 
temperature. 
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